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Abstract 

Questions of flux regulation in biological cells raise renewed interest in the narrow escape 
problem. The often inadequate expansions of the narrow escape time are due to a not so 
well known fact that the boundary singularity of Green's function for Poisson's equation with 
Neumann and mixed Dirichlet-Neumann boundary conditions in three-dimensions contains a 
logarithmic singularity. Using this fact, we find the second term in the expansion of the narrow 
escape time and in the expansion of the principal eigenvalue of the Laplace equation with mixed 
Dirichlet-Neumann boundary conditions, with small Dirichlet and large Neumann parts. We 
also find the leakage flux of Brownian particles that diffuse from a source to an absorbing target 
on a reflecting boundary of a domain, if a small perforation is made in the reflecting boundary. 

PACS numbers: 05.40.- Jc, 87.10.-e 



*Electronic address: amits@math.princeton.edu 
^Electronic address: schuss @ p os t . t au. ac . il 
^Electronic address: david.holcman@weizmann.ac.il 



1 



I. INTRODUCTION 



The narrow escape problem in diffusion theory, which goes back to Lord Rayleigh 
is to calculate the mean first passage time, also called the narrow escape time (NET), of 
a Brownian particle to a small absorbing window on the otherwise reflecting boundary of 
a bounded domain. The renewed interest in the small hole problem is due to its relevance 
in molecular biology and biophysics. The small hole often represents a small target on a 
cellular membrane, such as a protein channel, which is a target for ions Q , a receptor for 
neurotransmitter molecules in a neuronal synapse [3] , a narrow neck in the neuronal spine, 
which is a target for calcium ions 4|, and so on. The physiological role of the small hole 
is often to regulate flux, which carries a physiological signal. For example, the NMDA 
channels in the post synaptic membrane in the neuronal cleft are small targets for diffusing 
glutamate molecules released from a vesicle at the pre synaptic membrane. The leakage 
problem here is to find the fraction of the released molecules that reach the channels before 
being irreversibly absorbed by the surrounding medium (e.g., glia transporters) j^], {f| 



(see also http://en.wikipedia.org/wiki/Chemical_synapse). The position and the number 
of the NMDA and AMPA receptors regulate synaptic transmission and is believed to be 
a part of coding memory 

The narrow escape problem is connected to that of calculating the principal eigenvalue 
of the mixed Dirichlet-Neumann problem for the Laplace equation in a bounded domain, 
whose Dirichlet boundary is only a small patch on the otherwise Neumann boundary. 
Specifically, the principal eigenvalue is asymptotically the reciprocal of the narrow escape 
time in the limit of shrinking patch. 

The recent history of the problem begins with the work of Ward, Keller, Henshaw, Van 



De Velde, Kolokolnikov, and Titworth 
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ll| on the principal eigenvalue and is based 
on boundary layer theory and matched asymptotics, in which the boundary layer equation 



y 



jy Weber in 1873 
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3|. The 



on the 



is the classical electrified disk problem, solved explicit 

n 

work of Holcman, Singer, Schuss, and Eisenberg [14|, 
NET for diffusion with and without a force field and for several small windows and its 
applications in biology, is based on the known structure of the singularity of Neumann's 



2 



function at the boundary 



13 



22 



23 



241 ] and on the Helmholtz integral equation 25] 



(see [26(1). The most recent work of Benichou and Voituriez [27( on the NET in diffusion 
and anomalous diffusion finds the dependence of the NET on the initial point inside the 
boundary layer and finds the scaling laws for sub- diffusions. In these papers the leading 
term in the asymptotic expansion was calculated in the shrinking window limit. 

Neither the second term, nor its order of magnitude were calculated for the three 
dimensional problem, except in the case of a spherical domain with a small circular ab- 
sorbing window, where an explicit solution was constructed by a generalization of Collins' 



method (an error in the coefficient of the second term, given in jl5|, is corrected here). 
The difficulty in finding, or even estimating, the second term can be attributed to the 
practically unknown (to mathematicians and physicists) structure of the singularity of 



Neumann's function on the boundary. While classica" 



tions and in classical mathematical physics 
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22 



23 



texts in partial differential equa- 



24j mention only the leading order 



singularity of the Newtonian potential and a regular correction, 23] shows (in an exercise) 
that Neumann's function for a sphere has a logarithmic singularity at the boundary. The 
logarithmic boundary singularity of Neumann's function for the Laplace equation in a 
general regular domain seems to have been discovered by Popov 28] and elaborated by 
Silbergleit, Mandel, and Nemenman 29j] (which cites neither 23J] nor 281]) . 

Another small window problem is that of a leaky conductor of Brownian particles, 
which is a bounded domain with a source of particles on the boundary or in the interior, 
and a (big) target, which is an absorbing part of the boundary. The remaining boundary 
is reflecting. If the boundary has a small absorbing patch (a hole), some of the Brownian 
particles may leak out and never make it to the big absorbing target. The calculation 
of the leakage flux is not the same as that in the narrow escape problem, because the 
total flux on the boundary remains bounded as the small hole shrinks. The calculation 
of the leakage flux was attempted in 30] for diffusion in a flat cylinder with a source at 
the reflecting top and a small absorbing window at the reflecting bottom, and absorbing 
lateral envelop. The three-dimensional diffusion in the cylinder was assumed to be well 
approximated by radial diffusion in a circular disk. 



3 



In this paper, we find the structure of the boundary singularity of the Neumann func- 
tion for the Poisson equation and of the Green-Neumann function for the mixed problem 
(with Dirichlet and Neumann boundary conditions) in a general bounded domain Q, whose 



28 



29j . We find 



boundary dQ is sufficiently smooth. Our calculations use the method of 
that for z G dQ, y G Q U dfl, the structure of the Neumann function (in dimensionless 
variables) is 

1 1 



N(y,z) 



[L(z) + N(z)] In \y-z\+ v s (y, z), 



1) 



2iv\y — z\ 8n 

where L(z) and N(z) are the principal curvatures of dQ at z and vg(y,z) is a bounded 



function of x, y in Q. If Q is a ball of radius R, the above mentioned result of Kellog 23] 

is recovered, because L(z) = N(z) = —. 

R 

We find that the NET through a circular disk of (dimensionless) radius a, centered at 
on the boundary, is 



Et 





Q 




AaD 


[ L(0)+N(0) 
2tt 


- a log a + o(a log a) 



(2) 



where D is the diffusion coefficient. If Q is a ball of radius R, then 

101 



Et 



AaD 



a , R 
1 + — log — + o 
irR a 



The result ([3]) corrects that given in 
straightforward manner, as in 151 ]. 



(3) 



15]. The case of an elliptic window is handled in a 



The principal eigenvalue Ai(a) ~ — — of the Laplace equation in Q with Dirichlet 

Et 

conditions given on a circular disk of dimensionless radius a and Neumann boundary 
conditions elsewhere has the asymptotic expansion 
AaD 



Ai(a 

The result flU 
was not given in 



' L(0)+N(0) , . 1 , 

1 H a log a + o{a log a) 

2tc 



for a — > 0. 



(4) 



provides the missing second term and estimate of the remainder, which 

iy,0,0,[i5]. 



For a leaky conductor, we find that the leakage flux through a circular hole of small 
(dimensionless) radius e, centered at 0, is 



J £ = 4eDu (0) (1 + O(eloge)) 
4 



(5) 



where u (0) is the solution of the reduced problem (without the leak) at the hole. 

Equation (jSJ) can be viewed as a generalization of (jlj) in the sense that the factor \ 
in (jlj) can be interpreted as the uniform concentration of the Brownian particle in Q. The 
uniform concentration is the solution of the stationary diffusion equation problem with 
Neumann conditions on the entire boundary, which is the reduced problem for narrow 
escape. Thus the concentration Uq(x) is a generalization of the fixed concentration 
in flU). 



II. THE SINGULARITY OF NEUMANN'S FUNCTION 

Consider a bounded domain f2 C M 3 , given by f2 = {x G 1R 3 : F(x) < 0}, where 
F G C 3 (R 3 ). Our purpose is to determine the singularity of Green's function for the 
Laplace equation in Q with Neumann boundary conditions (called Neumann's function) 
and of Green's function for the mixed Dirichlet and Neumann boundary conditions. 

The Neumann function N(x, y) for this domain is the solution of the boundary value 
problem 

A x N(x,y) = -5(x-y) + -^-, for x,yeQ (6) 



dN(x,y) 



dv. 



0, for xedQ, yeQ, (7) 



x 



where v{x) is the outer unit normal to the boundary dQ. If x or y (or both) are in 

dQ, then only a half of any sufficiently small ball about a boundary point is contained 

in Q, which means that the singularity of Neumann's function is — ; r. Therefore 

2n\x — y\ 

Neumann's function for y G dQ is written as 

N(x,y) = — — : + v{x,y), (8) 

27T\x — y\ 



where v (x, y) satisfies 



A x v{x,y) = — for x e Q, y e dQ (9) 



5 



and the boundary condition 

dv(x,y) 1 u(x) ■ (x — y) 



for x,y G dVl. 



(10) 



dvx 27r I a; — y\ 3 
Green's identity requires the evaluation of two integrals. The first is the volume integral 
which by ([6]) is 

/ [N(x, y) A x v(x, z) - v(x, z) A x N(x, y)} dx = / N(x, y)—- dx + v(y, z) 



f v(x,z)dx, 



and the second is the surface integral, which by (JTj) is 



<)<1 



dv(x,z) dN(x,y) 
N(x,y) — v(x,z) 



8V; 



x 



x 



dS; 



X 



1 



2n\x — y\ 



+ v{x,y) 



v\x) ■ [X — z) 
2tt\x — z\ 3 



dS X - 



Thus, for z G dQ Green's identity gives 

v(y,z)= — j^j- J [N(x, y) - v(x, z)\ dx 



1 



+ v(x,y) 



V\X) ■ [X — z) 

2rr\x — z\ 3 



dS x . 



11) 



an L 27r F - V\ 

To determine the singularity of this integral when y approaches z, we use the method 
of successive approximations to expand v(x, y) as 



v(x, y) ~ v (x, y) + v^x, y) + v 2 (x, y) + . . . 



12) 



where Vi + i(x, y) is more regular than Vi(x, y) (see 29|). For y or z (or both) in dQ, the 
first term is the most singular part 

1 f v{x) ■ (x — z) 



v {y,z) 



4tt 2 J dn \x -y\\x- z\ 



■dS; 



X 



(13) 



To extract its dominant part, we reproduce here, for completeness, the analysis of 28 j 
with only minor modifications. We consider z G dQ and assume that the boundary near 



6 



z is sufficiently smooth. Moving the origin to z, we set z = 0. Taking a sufficiently small 
patch dflz about z, we assume that it can be projected orthogonally onto a circular disk 
D a of radius a in the tangent plane to dVL at z. We can assume, therefore, that dVLz can 
be represented as 

£3 = fz(xi,x 2 ) = -L(z)xl + -N(z)xl + o(xl + xl) for (x 1 ,x 2 )eD a . (14) 

If a is sufficiently small, then o(xf + x 2 ) <C Lx^ + iVx 2 ,. This canonical representation 
(|14p assumes that dflz has at least one non-zero curvature and that the quadratic part 
in Taylor's expansion of f(xi,x 2 ) about the origin is represented in principal axes. 

The asymptotically dominant part as y — > z is determined by the integral over the 
patch dVtz-i which we write as 

/ x l _ f »( x ) -xdS x , 15 . 

V ° Vl ~ 47r2 Jaa z y/( Xl - yi f + {x 2 - y 2 ) 2 + (x 3 - y 3 ) 2 [x\ + x 2 2 + x 2 f 2 

In the representation (|14p 

(L(z)x u N(z)x 2 , -1) + o(^/x\ + x 2 2 ) 



u{x) 



V\X) ■ X 



y/1 + L 2 (z)x 2 + N 2 (z)x 2 

L{z)x\ + N(z)x\ — x 3 
/l + L 2 (z)x 2 + N 2 {z)xj 



dS x = y/1 + |V/ z | 2 dxi dx 2 ~ y 1 + L 2 (z)x 2 + A^ 2 (z)x^ cfci cfe 2 , 

so that 

v ( 0) ~ — / ( L ( z ) x i + -^"(z)^ 2 . ~ ^3) dxi dx 2 

W ° y ' ~ 47r2 /d» v /(a; 1 -y 1 ) 2 + (x 2 -y 2 ) 2 + (x3-y3) 2 [xf + x 2 + a; 2 ] 372 ' 

The patch is represented in polar coordinates in D a as 

jf L ( z ) , . , iV(z)^ 2 



(xi,X2,x 3 ) = cos 0, r sin 0, r f - cos </H — sin + o(l) J ), (17) 

so transforming y into spherical coordinates 

(2/1,1/2,2/3) = |y|(sin6»cos0 o , sinflsin^o, COS6 1 ). 



we can write (TT6"j) as 



v (y,O) ~ i(\y\A,o) 



where 

[^Lr 2 cos 2 + \Nr 2 sin 2 + o(r 2 )] r dr 



n\y\ 



1 /2 

o [r 2 + \y\ 2 — 2r\y\ sin cos (0 — <f> ) + 0(r 2 |^/| + r 4 )] 7 [r 2 + 0(r 



4\l3/2 



~ [Lcos 2 + iVsin 2 0] / m . (19) 

2 Jo [r 2 + liyl 2 — 2r\v\ sin 9 cos(0 — 0n)l 



'o [r 2 + |y | 2 — 2r\y\ sin # cos(0 — 0o)] 
Integration with respect to r gives 



o + \y\ 2 — 2r\y\ sin 6* cos(0 — 0o)] 



1/2 



a — 


y 


sin # cos(0 - 


- < 


i>o) + V a2 + 


y 


2 -2a 


y 


sin 9 cos(0 — 0o ) 




y 


(1 — sin# cos(<^ 


>-0 O )) 





In 



ln^- + 0(l), 
I V I 



for y 7^ 0. It follows from (fT8j) that for y ^ z the leading order singularity is 

u (y, z) ~ ^-{L{z) + iV(z)) In + 0(1). (20) 

o7r \y — z\ 



For further analysis of the 0(1) term, see [29]. 

The canonical representation ( JT4l) of a hemisphere of (dimensionless) radius R at the 

south pole is x 3 = R — \JR 2 — (x\ + so = iV(z) = — . Therefore, for \z\ = R, 

R 

N(y, z) = + J- In + 0(1), (21) 

2-K\y — z\ AttR \y — z\ 



in agreement with 



23 



p. 247, Exercise 4]. 



III. APPLICATION TO THE NARROW ESCAPE PROBLEM 



A. Escape through a small circular hole 



As mentioned in the Introduction, the narrow escape problem 



16 



13, 



flfl 



10 



n 



14 



15 



181 ] is to calculate the mean escape time of a Brownian particle from a bounded 



domain Q, whose boundary is reflecting, except for a small absorbing patch (or patches 
0, 0) 9Q a . We assume here that dVt a is a circular disk of radius a <C l^] 1 / 3 and that a 
ball of radius R^> a can be rolled on dQ inside Q. This means that there are no narrow 
passages in f2. We denote <9f2 r = <9f2 — dVl a and e = a/l^] 1 / 3 and investigate the limit 
e — > 0. We assume that all coordinates have been scaled with |fi| 1//3 , so that all variables 
and parameters are dimensionless. 

The MFPT u(x) from a point x G O to dQ a is the solution of the mixed boundary 
value problem 

Au(x) = -—, for x G Q (22) 
u(x) = for xe dtt a (23) 
= for xe dQ r , (24) 

du x 

where D is the diffusion coefficient. The compatibility condition, 

is obtained by integrating ( [221 over Q and using (1231 and ([241 . 

Green's identity and the boundary conditions (JTj), ( 1231) . and ( 1241 give 

/ ^(*,y)daj= / N(x,y)?^dS x + C, (26) 



where 



C = / u(x)<ix. (27) 



Following the argument in 



151 ]. we note that iV(«c, y) is an integrable function independent 



of dQ a , whose integral is uniformly bounded, whereas C — > oo as s — > 0. Setting g(x) = 
dii(x\ 

— for x G dQ a and using the boundary condition (|23|) . we obtain from (|26|) the 

integral equation for the flux density g(x) in dQ a , 

N(x,y)g(x)dS x = -C for y G dVt a , (28) 



which, in view of (JSJ), (1201) now becomes the generalized Helmholtz equation [25|, []J 
1 



an a 



2n\x — y\ 



+ H(x,y) log 1^-1/1+0(1; 



dS x = -C for ye9fl fll 



(29) 



1 



H(x,y) = -—[L(y)+N(y))~-— 
Sir o7T 



[L(0) + iV(0)], for x,y & dQ a for e -> 0, 



where L(0),iV(0) are the principal curvatures at the center of dfl a - To solve (|29|) . we 
expand p(a;) = go (a?) + 9i( x ) + 02 (a?) H — ■ > where gr i+ i(aj) <C <7i(a;) for e — > and choose 

-2C 



9o(x) 



air\ 1 



Hi 



(30) 



It was shown in |l|], [26|, 15] that if dfl a is a circular disk of radius a, then 



0o (a>) 



27r ./an, I* - !/ 
It follows that <7i(a?) satisfies the integral equation 



dS x = C for all y e dVt a . 



(31) 



1 



2 ?r Jm a x-y 



9i(x) dS = 2C_ f H(x,y)\og\x-y\ 



an 



an a 



dS x . 



\x\ 



(32) 



Setting y = ar], x = a£, and changing to polar coordinates in the integral on the right 
hand side of (1321 . we obtain 



1 



2tt J ma \x-y 



g x {x) 2Ca 2 
aox 



2 /■27T 



an 



dO 



1 H(a£, at]) [log a + log |£ - 



r dr, (33) 



which gives in the limit e — > (e.g., keeping fixed and a — > 0) that 



dS'a; = — a log a + o(a log a). 



2tt Jan a |aJ - y 



2tt 



(34) 



As in the pair f[3"01 . f[3"Tj) . we obtain that 

-C[L(0)+iV(0)] 



9i{x 



7T 2 \ 1 



log a + o(loga). 



(35) 



\x\ 
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Finally, to determine the asymptotic value of the constant C, we recall that g(x) 
and use in ff25l) the approximation 

g(x) ~ g (x) +gi(x) ~ - 



du(x) 
du 



x 



-2C 



\x\ 



\ , L(0) + N(0) 
1 H a log a 

Z7T 



(36) 



V a z 

We obtain the narrow escape time Et = C (in dimensionless variables) as 
Et 










AaD 


r L(o) + iv(o) 

2?r 


- a log a + o(a log a) 



(37) 



The principal eigenvalue Ai(a) ~ — — of the Laplace equation in Q with the mixed 

Et 



Dirichlet-Neumann boundary conditions 

AaD 



Ai(a) 

The result fl38l) provides t 



, L(0) + iV(0) . . . 

1 H a log a + o(a log a 

27T 



has the asymptotic expansion for e — > 

(38) 



was not given in 



13,0, 



ie missing second term and estimate of the remainder, which 



11 



15|- 



If Q is a ball of radius R, then L(0) + N(0) = — and the narrow escape time Et = C 

R 



is given (in dimensional variables) by 



Et 



mi 



AaD 



i a i R 
1 log - 

ttR a 









-)1 




a /_ 



AaD 



R 



1 H log ho — log — 



nR 



R 



R 



(.39) 



The result f l39|) corrects that given in [15|]. Specifically, equation (3.52) in [lSj is missing 
the factor 1/tt of equation ( 1391) . which should have been carried from eq.(3.51) in 15] . 
The case of an elliptic window is handled in a straightforward manner, as in [151 ]. 



B. Leakage in a conductor of Brownian particles 

A conductor of Brownian particles is a bounded domain Q, with a source of particles 
on the boundary or in the interior and a target, which is an absorbing part dQ a of dQ. 
The remaining boundary dQ r is reflecting. Some of the Brownian particles may leak out 
of Q if dQ r contains a small absorbing hole S(e). The calculation of the leakage flux is not 

11 



the same as that in the narrow escape problem, because the total flux on the boundary 
remains bounded as the small hole shrinks. Our purpose is to find the portion that leaks 
through the small hole out of the total flux. 

The (dimensionless) stationary density u(x) of the Brownian particles satisfies the 
mixed boundary value problem 

DAu(x) = for xett 



-D 



du(x) 
dv 



du(x) 
dp 



(40) 



\x 



u ( x )\ 9 n a = u ( x )\s(s) = °> 



where <f)(x) is the flux density of the source on the boundary. Next, we derive an asymp 
totic expression for the flux through S(e), 

du(x 



Jr = D 



S(e) 



dv 



dSx, 



(41) 



in terms of the solution uq(x) of the reduced problem (without S(e)), thus avoiding the 
need to construct boundary layers. First, we find the flux of each eigenfunction and 
then, using eigenfunction expansion, we calculate J £ . Every eigenfunction u s (x) of the 
homogeneous problem (140!) satisfies 

-DAu £ (x) = X(e)u £ (x) for xett (42) 



du £ {x) 
dp 



for x e dVl, U dtt T 



(43) 



u £ {x) = for x G S(e) U dVl a 
The matched asymptotics method of 



(44) 



ll| gives the expansion of the eigenvalues 
A(e) = A(0) + Aie + o(e), (45) 
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where A(0) is the eigenvalue of the reduced problem (for fl without any small holes). 

We define the reduced Green function (without the small hole) as the solution of the 
mixed boundary value problem with D — 1, 

-AG(x,y) = 8{x-y) for x,yeQ (46) 
dG 

—(x,y) = for x e dtt s U Q r , y G ft (47) 
av 

G(x, y) = 0, for x G dtt a , yen. (48) 

Multiplying (146]) by u £ (y) and integrating over Q, we get 

^( £ ) f ni \ I \ -J f ru s9u e (y 



u £ (x) = -y- / G(x,y)u £ (y)dy + / G(x,y)—^dSy. (49) 



In view of the boundary condition (jHJ), we get from (j4"9~l) for all x G S{e) 

'S(e) 

The integral on the left hand side of (!50|) can be expanded about the center of 5(e) in 
the form 



/ G{x, y)u e (y) dy = - f G{x, y )^MldSy. (50) 

Jn JS(e) av 



[ \(e)G(x,y)u e {y)dy = G (e) + O(\x\) for x G 3{e), (51) 
Jn 

where the origin is assumed to be in the center of S(e) and the (21,22) plane is that of 
S(e). 

As in Section IIII Al Green's function for the mixed boundary value problem has the 
form 

G 0> y) = ttt- — 1 + H ( x > y) lo § \ x - v\ + v s( x > i/)> ( 52 ) 

2ty\x — y\ 

for x G dfl, y G QUdQ, where H(x, y) depends locally on the curvatures of the boundary 
and vs(x, y) is a continuous function of x, y G fl and on dfl. We assume that H(x, y) is 
bounded. Using fl52l) and the expansion (|5ip in (1501) . we obtain the Helmholtz equation 



^ + o(|*|) 



5(e) 



+ #(cc, y) log I a: - y | + u s (aj, y) 



27r|a; — y 



—^-dSy.(53) 
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The leading order singularity of G(x,y) and (1511) suggest the expansion 
du £ (y) C (e) 



<9z/ 



I V I 



+ 0(12/1) for ye S(e) 



where Co(e) is yet an undetermined coefficient, that is, 



G {e) 
D 



+ 0(\x\) 



S(e) 



2tt\x — y\ 



H(x, y) log \x-y\ + v s (x, y) 



C (e) 



f - 



\v 



+ 0(\y\) 



dS. 



y 



which reduces at x = to 

GJe) -C a (e)TT£ 



D 



+ I 0(\y\) 

S(e) 



2n\y\ 



+ H(0,y) \og\y\ 



dSy — 



C Q (e) [H(0,y)log\y\+v s (0,y)} dS y 



It follows that 



so that 



5(e) 



G,(e) 



\y\ 



+ [ 0(\y\)dS y 

JS(e) 



D 



C (e 



(y + 0(£ 2 loge)) C (e) + O(e 2 logs) 
G (e) + O(e 2 \oge) 



D 



tie 



+ 0(e 2 loge) 



Now, f[5"4"|) gives the flux through S'(e) as 



G? (e) + 0(e 2 \oge) f dS- 



0(e 2 \oge) Js(e) L_\y_ 



+ 



4s 



0(\y\)dS y 

S(e) 

G (e) + O(e 2 \oge) 



1 + 0(e\oge) 
14 



+ 0{e 2 \oge). 



To determine G (e), we integrate (135]) . to get the total flux condition 

X(e) [ u £ (x) dx = D I dS y + D [ ^ dS y . 

Jq Js{e) v v " Jan a ov 

We also recall that (jl8|) implies that 

/ G(x,y)^^-dS y = for x G dQ a , 

hence, using equations floDl) and flo^]) . we get the two equations 

G (e) + O( £ 2 log£) 



A(e) / u e (a:)dx = 4e " ' v ^ + O e 2 log£ 

Vn l + 0(eloge) 



Jau a dp 



A(e) / G(0,y)u e (y)dy = G (e) 

This gives 

A(e) / w e (cc) da; = 
in 



4 £ A( £ ) / G(0,?/K(y)^ + O(£ 2 log£) 

J * + 0(sHo g e) + D 

l + 0(elog£) J ma dv 



dSy 



Solving for A(e), we find that 



X(e) 



D J 


f du ' iy) dS v + 0(^o g e) 




f G(p,y)ue(y) dy + 0(e 2 \oge) 



D 



an a 



du £ (y) j C 



( As [ G(0,y)u £ (y)dy 
1 + ^Ja 



uJx) dx 



\ 



+ 0(e 2 \oge). 



u £ (x) dx 
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Note that 

f du £ (y) 

./,■«„ = A(0) + O(e), (61) 

u s (x) dx 

n 

due to the contribution of the boundary layer near S(e). 

Obviously, u £ — > tt as £ — ► 0, where u is the corresponding eigenfunction of the 
reduced problem (in the absence of the small hole, see also [sj]), so 



lim / G(x,y)u £ (y)dy = [ G{x,y)u (y) dy, lim f H^lJL dSy = [ dSy. 
£ ^°Jn Jn e ^°Jdn a ov J &n a av 

Therefore, using (159|) - (|6T]) in (1571) . we find that the flux of u e (x) through the small hole is 

J{e) = -D [ ?^ldS y = Ae\(0) [ G(0, y)u (y) dy + 0(e 2 \oge) 
JS(e) av ' Jn 

= 4:6 Du (0) + 0(e 2 logs). (62) 

Finally, expanding the solution u(x) of (j40p in eigenfunctions, we obtain from fl62l) 

J £ = Ae Du (0)(1 + O(e logs)), (63) 

where uq(x) is the solution of the reduced problem ()40l) . In dimensional variables, we 
obtain 

J £ = 4aD Po (0) + O ( log , (64) 

where po(0) is the value of the reduced stationary density (without the perforation) at 
the hole. 



IV. SUMMARY AND DISCUSSION 

The main results of this paper are (i) the explicit calculation of the second term in the 
expansion of the NET, which can be quite significant, and which also provides a bound for 
the remainder in the expansion; (ii) an explicit expression for the leakage flux through a 
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small opening in the impermeable envelope of a conductor of ions. The leakage is often a 
key control mechanism of physiological function, such as in the synaptic cleft of a neuron, 
as mentioned in the Introduction. The leakage formula (|64p can give explicit expressions 
for the flux when the reduced problem is explicitly solvable, e.g., in simple geometries. If 
there are several leaks, at Xi, then ( |64l) gives 

J £ = 4aD ]T po(^) + O (r^j- 3 log ° ) , (65) 

which demonstrates the role of clustering or un-clustering of the leaks in regulating flux 
[2^. Specific applications of the results of this paper to molecular biology and 
biophysics will be published in a separate paper. 
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